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Atmospherics are assumed to be spike shaped with exponential upsurge and decay; 
their occurrence in time is Poisson-distributed. The distribution of the peak value is shown 
to be lognormal. Under these conditions, the amplitude distribution of atmospheric noise 
is derived; it may be well approximated by the distribution of the sum of a Rayleigh vector 
and a vector with lognormally distributed amplitude. The resulting distribution is uniquely 
determined by three physically meaningful parameters and is in very good agreement with 
measured distributions. 

1. Introduction 

The probability distribution of the amplitude of atmospheric radio noise is an important 
basic factor for the design of radio systems; it may also be used as an effective instrument for 
geophysical research of atmospheric discharges by which this noise is caused. The statistical 
distribution of the amplitude of atmospheric radio noise is therefore being measured in observa- 
tories all over the world; so far it has, however, not been possible to give a thorough and 
practically useful theoretical explanation of this distribution. 

It has been found empirically that the distribution 

P(E>R) (1.1) 

where E is the instantaneous value of the envelope of the HF atmospheric noise and R is a 
selected level of voltage or field strength at the receiving point, approaches a Rayleigh dis- 
tribution for small amplitudes (high probabilities). This is evidently caused by the interference 
of many mutually overlapping small atmospherics whose phase is uniformly distributed over a 
basic phase cycle, for these atmospherics are emitted by independent sources at distances of 
many wavelengths from the receiver and from each other. For large amplitudes (low proba- 
bilities) the distribution deviates considerably from a Rayleigh distribution and may be well 
represented by a lognormal distribution [Horner and Harwood, 1956; Foldes, 19(30]. This 
part of the distribution curve is evidently due to relatively rare, strong atmospherics that do 
not overlap in time. 

It has also been found empirically [Crichlow et al., 1960] that when the experimentally 
measured distribution (1.1) is mapped on Rayleigh paper (cf, appendix), the lower portion of 
the curve (small R) is a straight line with slope corresponding to the Rayleigh distribution and 
the upper portion of the curve (large E, small P) may be adequately represented by a straight 
line of larger slope than that corresponding to a Rayleigh distribution. As shown in the 
appendix (and fig. 1 ), a lognormal distribution will map on Rayleigh paper in this way. Between 
the two straight lines there is a transition that has been approximated by Crichlow et al. 
[1960] as a circular arc. 

The individual atmospherics that protrude from the permanent noise level are spike 
shaped and decay rapidly, exceeding the permanent noise level only for some milliseconds 
[Horner and Harwood, 1956]; this order of time is in agreement with the characteristics of 
lightning discharges [Watt and Maxwell, 1957]. The decay is evidently exponential and has 
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Figure 1. Log-normal distributions plotted on 
Rayleigh paper. 

a small time constant. The buildup is also exponential with a still shorter time constant. 
The purpose of this paper is to derive the theoretical probability distribution (1.1). In 
section 2 we shall first derive the probability distribution of the peak values of the individual 
atmospherics. In section 3 we derive some general probabilistic relations, in particular the 
distribution of the time interval separating a peak from the Mi following or the kth. preceding- 
peak assuming the number of peaks in a given time interval to be Poisson distributed. In 
section 4 we derive the amplitude distribution of the decaying Hh. preceding atmospheric and 
the rising kih future atmospheric; the time constants of rise and decay may be chosen indi- 
vidually and arbitrarily, so that we are working with a very general spike-shaped atmospheric. 
In section 5 the rms value of an individual atmospheric and of the entire atmospheric noise is 
derived and the required distribution is then found by summing over all k; the resulting distri- 
bution is the sum of a Rayleigh vector and a vector with lognormally distributed amplitude. 
In section 6 this distribution is compared with the experimentally measured distribution; in 
spite of some coarse approximations, the agreement is surprisingly good. As in the case of 
the empirically graphical method of Crichlow et al. [1960], the resulting distribution is uniquely 
given by three parameters; as one of these depends exclusively on the properties of atmospheric 
discharges and one predominantly on propagation conditions (the third drops out in normalizing 
to the rms value), the theoretical distribution derived in this paper should prove useful for 
geophysical research. 

2. Distribution of the Peak Value of Atmospherics 

The peak value E p of an atmospheric is determined by the peak power P of the atmos- 
pheric discharge emitting it, the distance of the discharge from the receiving point, and the 
attenuation A in the propagation of the atmospheric (all at the considered frequency) ; therefore 



^=^exp(-gM;) 



(2.1) 



where if is a constant determined by the system of units employed; the attenuation A is given 
by the sum of attenuations on the individual sections dj of the propagation path (whose total 
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length is d); dj is the attenuation constant for the jth section. We may thus write (2.1) as 

E p =e A =exp (-S 8j(/j+l In P -ln S d,+ln k) (2.2) 

which defines A. 

The first three terms of the exponent are random; we may assume that the fluctuations 
(the variance) of P and d will, on taking the logarithm of these quantities, be so much diminished 
that they will be small compared with the fluctuations (the variance) of the sum of attenuations 
as given by the first term of the exponent, which will thus predominate in determining the 
distribution of A. We further assume that the total propagation path d can be split up into a 
large number of sections dj on which the attenuations 8jdj may be considered independent with 
no section making a predominant contribution to the sum of attenuations. This assumption 
seems reasonable since most atmospherics are propagated from long distances via the ionosphere. 
Then A will be distributed normally, i.e., its probability density will be 

*^Hr«p r-^1 (2 - 3) 

where /x is the mean value of A and a 2 its variance. 
From (2.2) we have 

A=ln E p (2.4) 

which together with (2.3) yields the required distribution of E p : 



WP v (In E n — ju)H fo K v 

\{E P ) = — — ^=exp ^ ( 2 - 5 ) 



I 

<jE p ^2t 



and hence 



P(^>5)=A(fl)=J\(^)dE;=i[l-erf (k*_i»YJ (2.6) 



The distribution (2.5) or (2.6) is the logarithmic-normal distribution, which has been found 
experimentally to apply to strong atmospherics [Horner and Harwood, 1956; Foldes, 1960] but 
was not explained theoretically. The mean and mean square values of E v are from (2.5) 

<^) = exp( M +~) (2.7) 

(E/) = exp2( M +<r*). (2.8) 

From (2.8) and (2.6) we then have 

If we wish to express 11 in decibels, we obtain after elementary manipulations 

where 

A=20 log w 6 = 8.686. (2.1 I) 

Curves of (2.10) for various values of a are plotted on Rayleigh paper in figure I. 
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3. Probabilistic Relations 

In order not to interrupt our development later, we now derive two general relations that 
we shall need. 

Since atmospherics are emitted by mutually independent sources and the reception of an 
atmospheric is equally probable at any time, their number in a given time interval will be 
Poisson distributed, i.e., the probability of receiving n atmospherics in a time interval of 
length r is 

P{n] r)= ^r e ~ NT ' (3 - 1} 

where N is the mean number of atmospherics received per unit time. 

To find the probability density of the time interval T k separating a given peak from the 
kt\\ preceding peak, i.e., 

Vk(k)=Pr{t k < T k <t k +dt k ) (3.2) 

we note that this probability equals that of receiving k— 1 atmospherics during the time t k 
and one atmospheric during dt k ; since the two events are independent, we have 

p k (t k )=P(k-l; t k )P(l; cU k ) (3.3) 

or substituting from (3.1), 

AW -1 

ft^ = (E=i)] e ~ Ntk < 3 - 4 ) 

which is the first relation we shall require. 

The second concerns a random signal (vector) with independent amplitude and phase 

E=Ee**, (3.5) 

where we have suppressed the time factor exp (iut). Let the phase be distributed uniformly 
over an interval of length 2x, or more generally let the probability density of the phase, p(<f>), 
be such that for all in the interval 0<</><2tt 

±J^+2k^ = ^+e^); |€(*)|«^: (3.6) 

so that #>(<£) is equivalent to a uniform distribution in the interval (0, 2ir). A vector of the 
kind (3.5), whose phase distribution satisfies (3.6), is met so often in wave propagation that it 
deserves a name; because of its Uniformly Distributed Phase we shall call it a UDP vector. 
(A Rayleigh vector is the special case of a UDP vector with Rayleigh-distributed amplitude.) 
The sum of any number of UDP vectors 

A n a n 

U=Ue^=^E j =^E j e i ^ (3.7) 

is obviously again a UDP vector. Consider now the mean square value of a UDP- vector sum 

/ n n \ n 

=1 M£ 

But for J9^k, 

( e u<i> r <t> k ))=J- J^ J T e i( ^-^d<t>jdct> k =0. (3.9) 

Substituting (3.7) and (3.9) in (3.8) we find 






n a 

}=1 



)=jt (E 2 j\ (3.10) 
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regardless of the distributions (possibly all different) of the amplitudes E jm Thus in the 
interference of waves that may be represented by UDP vectors, the mean total power equals 
the sum of mean powers of each wave regardless of the amplitude distributions of the latter. 

The phase of an atmospheric is, among other factors, determined by the distance of the 
atmospheric discharge from the receiver. Since this distance is random and varies over a 
range much larger than one wavelength, (3.6) is satisfied and any atmospheric may be repre- 
sented by a UDP vector. 

4. Amplitude Distribution of the Signal Due to the icth Decaying Atmospheric 

We assume the shape of the envelope (amplitude) of an individual atmospheric attaining 
its peak value E v at time t of the form 



u (t)=* 



Vp)o«p(-^)forC>to 

JE p ) exp(^yort<t 



(4.1) 



i.e., in the shape of a spike as shown in figure 2a with arbitrary time constants of decay and 
upsurge a and b. This shape evidently corresponds to the real shape of atmospherics better 
than rectangular impulses (cf, fast film recordings by Horner and Harwood [1950]). 

However, in addition to this atmospheric u (t Q ) = (E p ) there will be present also an infinite 
number of decaying atmospherics u ly u 2 , . . . u k , . . ., that have attained their peak 
value (E p ) k in the past, t k time units ago (fig. 2b) ; the distribution of t k is given by (3.4). The 





Figure 2. Random sequence of atmospherics. 
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amplitude of the kth (past) atmospheric at time t is therefore 

(E p ) k exp (~)=g (4.2) 

A=(E P ) k , B=exp(^) (4.3) 



where 



are independent random variables. The distribution of A is X(^4) as given by (2.5). The 
distribution of B is found from (3.4) by transforming in accordance with (4.3): 

^(- B )^ ^_ 1 )Tga^V + l (1<5<oo). (4.4) 

The distribution of u is now found from (4.2) as the distribution of the ratio of two inde- 
pendent random variables [Gnedenko, 1954]: 



p(u*)= z\(uz)w(z)dz. 

Substituting In z=x we find 

/ v N k a k r k _ x r AT (In Uk -\-x- ^T\ i 
gW %(t-l)lcrVS Jo X ^ L"^ 2^~ J *' 

The mean square value of u k is 



(4.5) 



(4.6) 



/» CO 

(uj>= J u 2 k p(u k )du k (4.7) 

and may be found from (4.6) by interchanging the order of integration and the substitutions 
In u=y, (Na+2)x=t: 

{Uk} - (Na+2Y(k-l)\ J l e dL (4 - 8) 



The integral equals T(k) = (k— 1)1, so that we obtain 2 

2 _ A^Vexp[2((T 2 +)u)] 
^* (iVa+2) & 



(4.9) 



Besides the decaying atmospherics u k , we have to consider also the upsurging atmospherics 
s k attaining their peak value at a time t ~\-t k . For the upsurging atmospherics the calculation 
is of course the same except that the time constant a must be replaced by b so that 

jyy exp [2 (,*+„)] 

W (Nb+2)* { • } 

5. Amplitude Distribution of Atmospheric Noise 

The total signal U at time t is 

#=£ +ibVfib«* (5.1) 

/.- = ! fc=l 



2 More generally, we find by the above method 
of which (4.9) is a special case for ra=2. 
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where the circumflex accents again denote UDP vectors. Applying (3.10), we therefore find 
from (2.8), (4.9), and (4.10) 

<CP)=exp [2 (^+M)]{l+g (- N ~J+± (wM}- (5-2) 

Both series are geometric with sum Na/2 and Nb/2 respectively. Introducing the average 
time constant of upsurge and decay 

c= a ~¥ (5-3) 

we then have 

(U 2 ) = (l+Nc) exp [2 (<t 2 +m)]. (5.4) 

This is the mean square value of the field at a time t Q , the time of arrival (peak) of an 
atmospheric. The probability distribution of the amplitude at that lime is strongly dependenl 
on the value of Nc. It is easily shown from (5.1) and (5.4) that for 7VV>>1, which corresponds 
to the atmospherics overlapping in time at all but the highest levels, U is approximately Ray- 
leigh distributed; this contradicts the measured distribution. 

On the other hand, if iVc<^<0> the two series in (5.2) will be much smaller than unity; 
in this case the log-normally distributed quantity u will predominate and the overlapping 
decaying or upsurging atmospherics u k and s k will affect the resulting distribution only when 
u (and hence U) is small. This is in agreement with the measured distribution and we therefore 
henceforth assume 7Vc<C < 0. Also, we may safely assume 6<C<^a, i.e., that the time 
constant of the upsurge is much smaller than that of the decay (cf, recordings by Horner 
and Harwood [1956]). We thus assume 

Nc^~«l. (5.5) 

A 

It follows from (5.1) that U is a UDP vector composed of infinitely many UDP vectors, 
of which one (u ) has a mean square comparable to the mean square of the sum, whereas the 
mean squares of all the others are much smaller and these remaining vector terms will affect 
the distribution of i^only for small amplitudes of u {) and hence of U. The sum of these remain- 
ing vectors is that of a rapidly converging geometric series and thus the conditions for it to be 
a pure Rayleigh vector are not strictly satisfied [Beckmann, 1962a], though its distribution 
curve will asymptotically approach a Rayleigh distribution for small amplitudes and thus 
this sum may well be approximated by a Rayleigh vector. Thus U as in (5.1) will equal the 
sum of a lognormally distributed vector and a Rayleigh vector. The distribution of the 
amplitude of the sum is then given by [Beckmann, 1962b] 

W{U) = W e ~ U2 ' M j^ H^)e- ul/M I (2u U/M)dn (5.6) 

where Mis the mean square value of the Rayleigh vector and I is the modified Bessel function. 
It can be shown [Beckmann, 1962b] that for small values of U, (5.6) will approach a Rayleigh 
distribution, whilst for large values of U it will approach a lognormal distribution. 

However, it should be pointed out that (5.6) is the distribution of the amplitude at time 
t , i.e., at the time of arrival (peak) of an atmospheric. For any arbitrary time t (between two 
successive peaks) the amplitude is 

E=Ue~ t/a (5.7) 

(where the upsurge has been neglected since b<C<^^) and although U is random, once it has 
assumed a certain value, E is governed by the deterministic exponential law (5.7) between 

two consecutive atmospherics. Thus E at any time is also the sum of a lognormal and a 
Rayleigh vector (both are UDP vectors). If the time interval separating two successive 
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atmospherics is r, then from (5.7) 

(E 2 \ T )=(U 2 ) - C e- 2l/a dt=£- (U 2 ){l-e~ 2 ^) (5.8) 

T JO *T 

or using (5.5), 

(E 2 \r)=(U 2 ) C -(l-e-^). (5.9) 

Now r is itself a random quantity with probability density Ne~ Nr , i.e., (3.4) for k=l; 
applying the Theorem of Total Probability to (5.9) we therefore have 

(E 2 ) = {U 2 )Ncrr- 1 (l-e- r/c )e- Nr dr=(U 2 )Ncln^^ (5.10) 

or from (5.4) and (5.5) 

(E 2 ) = (l+Nc)Nce 2 ^ 2+ ^ In^^^Nce 2 ^ 2 ^ In (1/iVc). (5.11) 

Similarly, the mean square value of E t =Ui exp (—t/a) is 

(Ei) = (ui>Nchx(l/Nc), (5.12) 



in particular, 



(El)=Nce 2( °^ \n(l/Nc). (5.13) 



It follows from (5.11) and (5.12) that the predominant term in E is E , whilst the sum of 
all remaining terms 

A oo A 

i=\ 

which is Rayleigh distributed, has a mean square value of only 

(Z 2 ) = (E 2 ) - (El) =Nc(El) =N 2 c 2 e 2( ° 2 ^ In (l/Nc) . 

To find the distribution of E we use E =U e~ t/a where U is lognormal, t is uniform 
from to r and r is distributed according to (3.4) for k=l; hence 

(5.14) 

This integral cannot be evaluated in a closed form, but approximate calculations show 
that the distribution of E Q remains roughly lognormal, especially for large values of E (where 
the exponential decay distorts the distribution of U less than for low values attained after a 
longer time r). Thus if we replace the random variable r by its mean value ( T ) = l/N and 

use (5.5), we find 

-E-o \ ay l J 

which for large values of the argument (i.e., large E ) will asymptotically tend to 

,_,, 27W 1 r (In g -M)H 

*™ -ef; • ^E~m exp L — 2^-J 

thus retaining the character of a lognormal distribution. Alternatively, one may expand the 

square bracket in (5.14) in powers of r; this again yields a lognormal distribution that may be 

factored out of the integral. Other approximations of (5.14) again retain a lognormal character. 

We therefore approximate (5.14) by a lognormal distribution with mean square (5.13). 
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Since (5.13) may be written as 

E* =exp ^W'+^-hx Nc+\\n\n (1/Nc)}\, 
the required distribution may be taken as log-normal with parameters <r'=a, and 

M '= M +I In Nc+l lnln (1/iVc). 

The distribution of A 7 will therefore be the sum of this lognormally distributed vector with 
amplitude E and the Rayleigh-distributed vector Z. 

The required distribution of E is therefore in analogy with (5.6) given by 

P{E)=2 M e ~ E2IM f ~ Po( E o)^ E ° 6lM M^U E/M)dE (5.15) 

where M=(Z 2 ) and p (E ) = \(E ;a;fjL / ) t 

Since, in our case, El is much larger (1/iVc times) than M, (5.15) turns into p (E ) for 
large E and into a Rayleigh distribution for small E Q . Normalizing with respect to the root- 
mean-square value of E, which from (5.11) is 



#nns « jNc In (l/Nc)e oi+ll =(E ) mB9 (5.16) 

these two distributions are therefore 

for large R, and 

p(E>B)=e Wc (5.18) 

for small A*. 

For the intermediate (and all other) values of R we find by normalizing and integrating 

(5.15): 

T>( E ^li\ 2 C<IF [".IF FF-i ex-n T E *+ El Q"^o+OH f f 2J ^ ,r 1Q x 

p v^ >7; rivw^J/ I ° ° exp L~~iv<r -a?~J 7 <-w/ (5 - 19) 

This is the formal solution of our problem; the practically useful solution is given by the 
asymptotic distributions (5.17) and (5. IS). In comparing our results with the measured 
distribution (section 6) we have used (5.17) and (5.18). For the intermediate region, which 
turns out to be comparatively small, we have not gone to the trouble of computing the cumber- 
some integral (5.19), but have followed the example of Crichlow et al. [1960] in joining the 
curves of the two asymptotic distributions (5.17) and (5.18) by a circular arc (dashed section 
of curves in figures 3 to 5). 

6. Comparison With Experiment and Applications 

Crichlow et al. [1960] have made an empirical analysis of more than 100 amplitude- 
distribution curves of atmospheric noise measured in Colorado, and further curves measured 
in Alaska, Panama, Florida, and England. In all these cases they found that when the dis- 
tribution curve (1.1) is plotted on Rayleigh paper, the measured points lie for small R on a 
straight line with Rayleigh slope, which is in agreement with (5.18), and on a straight line with 
larger slope for large B; this is again in agreement with (5.17) for as shown in the appendix, 
the lognormal distribution will for large R plot as a straight line on Rayleigh paper and its 
slope is determined by the parameter a. 
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The position of the Rayleigh line is determined by a single point (its slope is already 
determined by the coordinate net). From (5.18) it follows that for 



we obtain P{EjE TTaa >R)-- 



= 37%. 3 



E=^Nc 



(6.1) 



Thus the parameter Nc determines the position of the Rayleigh line, a determines the 
slope of the high -voltage, low-probability line and the parameters o~, Nc, and jjl determine 
the rms value of the atmospheric noise E vms . From the above it follows that the distribution 
(5.17) to (5.19) must always be in agreement with the measured distributions analyzed by 
Crichlow et al. [I960]. Examples are shown in figures 3,4, and 5. The distribution agrees 
equally well with the results measured by Watt and Maxwell [1957]. 

The distribution derived here is uniquely determined by three parameters: Nc, a, and jjl. 
The first of these depends only on the properties of atmospheric discharges (N the number of 
discharges per unit time, a the time constant of decay) ; a is the standard deviation of the total 
attenuation, which in accordance with section 2 is almost exclusively determined by the proper- 
ties of the propagation path ("almost" because to some small extent it depends on the standard 
deviation of the logarithm of the peak power of the discharge) ; /z is the mean value of the total 
attenuation and is also predominantly determined by the properties of the propagation path 
(this parameter drops out in the normalization with respect to the rms value). 
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Figure 3. Amplitude-probability distribution of 
atmospheric radio noise. 
Circles: values measured by Crichlow et al., [1960] on 13.3 kc/s at 
Boulder, Colo., October 6, 1958. Full curve: distribution computed 
from (5.17) and (5.18) for o-=1.67, iVc=0.01. Broken curve inter- 
polated. 
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Figure 4. Amplitude-probability distribution of 
atmospheric radio noise. 
Circles: values measured by Crichlow et al., [1960] on 10 Mc/s at 
Boulder, Colo., September 10, 1958. Full curve: distribution 
computed from (5.17) and (5.18) for <r = 1.15, Nc=0.1. Broken curve 
interpolated. 



s This point is usually not marked on the P scale; a better policy is therefore to determine the point corresponding to P= 
as may be verified from (5.18), lies almost exactly 20 db below the rms value of the Kayleigh component (6.1). 
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Figure 5. Amplitude-probability distribution of 
atmospheric radio noise. 
Circles: values measured by Crichlow el al., [1960] on L13 kc/s at 
Boulder, Colo.. September 12, 1958. Full curve: distribution 
comp ui( m1 from (5.17) and (5.18) for <r = 1.74, Nc=0.008. Broken 
curve interpolated. 



Applications to geophysical research arc evident: a change of the above parameters will 
cause a change in the amplitude-distribution curves as measured at different receiving points, 
times, and under different conditions. Comparison of the measured curves with the theoretical 
curves should permit the assessment of the causes resulting in the particular measured 
distribution. 

Among other applications of the theoretical distribution, the following one seems hopeful: 
the variations of the level of atmospheric noise may be caused by a change in propagation 
conditions or lightning activity or both, and it has been an obstacle to geophysical research that 
the two causes are difficult to separate, The theoretical distribution derived above should 
enable the two causes to be separated: if there is an increase in lightning activity, the parameter 
A r c will increase, which will prolong the Rayleigh portion of the curve and raise it with respect 
to the rms value. If propagation conditions change, this will affect the parameter a and hence 
the slope of the high-voltage, low-probability portion of the curve. 

The problem of the amplitude distribution of radio signals scattered by ionized meteor 
trails is very similar to the problem treated here, since the occurrence of meteor trails, like 
that of atmospherics, forms a Poisson process in time and the echo (often) decays exponentially. 
The method used here may therefore in principle be applied to improve the theory of that 
distribution [Beckmann, 1962c]. The main difficulty is that the distribution of the peak 
ralues (dependent on the distribution of the mass of meteors) as used in astronomy lias no rms 
value. 

The method of calculating the overlap of time-dependent (decaying) phenomena whose 
occurrence is Poisson distributed in time as introduced in section 4 may have applications for 
queueing and inventory processes. 
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The main points of this paper were worked out at the Institute of Radio Engineering and 
Electronics, Czechoslovak Academy of Sciences. 4 It was rewritten in a revised form after 
A. D. Spaulding of the NBS Boulder laboratories had made a thorough study of the original 
version and pointed out some errors in it; the author wishes to express his sincere gratitude to 
A. D. Spaulding. The author is also grateful to Mrs. Charlotte Cranford, who typed the 
manuscript. 

7. Appendix. Mapping of Distribution Functions in Generalized 

Probability Coordinates 

Let the distribution of a random variable E be given by the relation 

P(E>R) = P(R). (A.l) 

We wish to map this distribution in Cartesian (£, rj) coordinates with scales ^(P), rioy(R) 
chosen in such a way that (A.l) will be mapped in these coordinates as a straight line. The 
moduls or scaling constants £ > Vo convert £ and rj to centimeters, inches, or other units of length 
and we set them equal to unity. 

Since (A.l) is always a nonincreasing function, the design of such coordinates is always 
possible. One of the two scales, usually r)(R), may be chosen arbitrarily (e.g., linear, logarith- 
mic, etc.), the other is then determined by the requirement 

V (R) = C+SZ(P) (A.2) 

where C and S are choosable constants; S is the slope of the required straight line. 

Let R(P) be the inverse function to P{R) as given by (A.l); then from (A.2) the required 
scale of the £ axis is 

«f)_2a£2» (A.3) 

which is the general solution of the problem. 
For the Rayleigh distribution 

PCE>i?)=exp(-g) (A.4) 

the inverse function is 

R(P)=E rm ^hi (1/P). (A.5) 

If we choose 

rj(R) =M log 10 R=M' In R (A.6) 

(where the constants M or M' may be absorbed into the modul rj 0y so that we may set either 
equal to unity), we obtain the scale £(P) of Rayleigh paper on which (A.4) will be mapped as a 
straight line with slope S, by substituting (A.5) in (A.6) and the result in (A.3): 

!CP)=2^1og[ln(l/P)] + a (A.7) 

where C x is a constant determining the position of the origin on the £ axis; we set it equal to 
zero 5 . The slope S for Rayleigh paper is usually chosen as S= — l or S=— %. The latter is 
used in figures 1 and 3 to 5 of this paper. 



4 URE-CSAV Inst. Rept. No. 26, 1962. 

5 Since either M or M' in (A.6) may be set equal to unity, the symbol "log" in (A.7) may mean either logio or In. 
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It may be shown by elementary methods that all functions of the kind 

P(B) = exp (-kin (A.8) 

and only these functions will map on Rayleigh paper (A. 6), (A. 7) as straight lines. The slope 
of such a line is 

sr / =— (A.9) 

m 

The curve that will result from mapping a general function P(R) on Rayleigh paper is 
found by substituting P(R) in (A.7) and eliminating R from (A. 7) and (A. 6), so that the curve 
V=v(0 w ill be obtained. For example, for the lognormal distribution normalized as in (2.9) 
we have 

P(5)=Jerfcs (A.10) 

where 

erfc x=l— erf x (A. 11) 

and 

In R+v 2 77 + a- 2 



7V2 (T-V2 



(A.12) 



In the second expression of (A.12) we have substituted from (A. 6) for M' = l. The re- 
lation between £ and 77 is therefore 

|=iln[ln2-lnerfc(^)} (A.13) 

This (or any other) function will be mapped on Rayleigh paper as a straight line in those 
regions where d^/drj=const. From (A.12) and (A.13) 

<^ = d! L dx = erf_ (A 14) 

drj dxdr\ Sa^j2Tr(\n 2— erfc x) erfc x 

This expression will be constant for x constant, i.e., from (A.12) for 

|lnP|«a 2 (A.15) 

making ^ ~o-/V2, which substituted in (A. 14) yields 1JS', where S' is the slope of this straight 
line. The ratio of S' to £ (where S is the slope for which the Rayleigh paper was designed), 
is therefore 

f =( 8 iy=°^» [in 2-ln erfc |] erfc J, (A.16) 

Thus the lognormal distribution will map on Rayleigh paper approximately as a straight 
line in the region where (A.15) holds, i.e., in the region of small probabilities for cr>l (cf, 
fig. 1, p. 724). 

One way of determining a from a measured amplitude-distribution of atmospheric noise 
is to plot the upper portion of the curve on Gaussian probability paper with the R axis given 
in decibels. If this part of the curve is lognormal, the points will lie on a straight line and o f 
may be determined, e.g., as the decibel difference between the 2.28 percent and 0.135 percent 
levels. The parameter a is then found from the relation 

o- / = 20alog 10 e = 8.686o-. 
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